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Abstract 



a . 

P ■ A Rayleigh-Benard cell has been designed to explore the Prandtl (Pr) dependence of tur- 

bulent convection in the cross-over range 0.7 < Pr < 21 and for the full range of soft 
and hard turbulences, up to Rayleigh number Ra ~ lO^"*^. The set-up benefits from the 
k>( I favourable characteristics of cryogenic helium-4 in fluid mechanics, in-situ fluid property 

■ measurements, and special care on thermometry and calorimetric instrumentation. The 

5^ \ cell is cylindrical with diameter / height = 0.5. The effective heat transfer Nu{Ra, Pr) has 

been measured with unprecedented accuracy for cryogenic turbulent convection experiments 
in this range of Rayleigh numbers. Spin-off of this study include improved fits of helium 
thermodynamics and viscosity properties. Three main results were found. First the Nu{Ra) 
dependence exhibits a bimodality of the flow with 4 — 7% difference in Nu for given Ra 
and Pr. Second, a systematic study of the side-wall influence reveals a measurable effect 
on the heat transfer. Third, the Nu{Pr) dependence is very small or null : the absolute 
value of the average logarithmic slope {dlnNu/ dlnPr)Ra is smaller than 0.03 in our range 
of Pr, which allows to disciminate between contradictory experiments [Ashkenazi et al, 
Phys. Rev. Lett. 83:3641 (1999)][Ahlers et al, Phys. Rev. Lett. 86:3320 (2001)]. 



1 INTRODUCTION 

Static equilibrium in a column of fluid corresponds to a balance between many parameters 
such as the weight, the pressure gradient, the temperature difference, etc. The occurrence 
of a small local perturbation can initiate a global convective motion. Rayleigh-Benard 
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convection is a reference configuration for convection studies: a fluid cylinder (height h, 
cross sectional area S) located between two horizontal plates is subjected to a temperature 
difference AT between the plates. In our case the upper plate is regulated at constant 
temperature and AT results from a constant heat flux Q applied at the bottom plate. For 
high enough AT the convective flow turns into a turbulent regime (for additional references, 
see for example ref. [1]). 

For a given fluid at mean temperature T with a mass density p, and for fixed geomet- 
rical conditions, the convective flow is characterized by one single parameter AT or the 
dimensionless Rayleigh number defined as: 

Ra=^^^ (1) 

UK 

In this expression: 

- a is the constant-pressure thermal expansion coefficient, 

- g is the gravity acceleration, 

- is the kinematic viscosity, k the thermal diffusivity ; their ratio is the Prandtl number: 

Pr = u/k (2) 
The Nusselt number gives the apparent thermal conductivity in the cell: 

where A is the fluid thermal conductivity. In a given cell, Nu should depend only of Ra 
and Pr. The influence of the adiabatic gradient on Ra and Nu have been compensated, 
with the exact correction formula[2] : 

ita — J^Oiuncorr AT 

AT - AT AT Q - XSATadia,/h 

^-t — adiab Q 

An additional correction due to the sidewall conductance is applied, according to the 
formula proposed in section 4.2. 

The unique properties of cryogenic ^He allow to control high Rayleigh numbers[3, 4, 5, 
6, 7]. In particular, these were used in Grenoble to reach Ra higher than 2.10^^ in "high" 
cells (aspect ratio 1/2, h = 20 cm) : in such conditions they observed for the first time 
the Kraichnan regime[6] (also called the ultimate regime) and its asymptotic limit[8], both 
predicted forty years ago [9]. 

Helium gives also the opportunity to easily vary the Prandtl number[10] over an unusual 
range. We have done a specific study of the Pr variation effect in a small size cell (called 
the mini-cell, aspect ratio 1/2, h — 2 cm), for 3.10^ < Ra < 10^^ corresponding to the 
soft and hard turbulence regimes where experimentally controlled Pr can be achieved 
independently of Ra (the lowest explored Ra in this work is below lO'' and the convection 
threshold was found to be the same as in the large cells, around 4.10^). This Ra range 
fully covers the turbulent regimes preceding the transition towards the ultimate regime and 
allows a comparison with other experiments using various fluids, and also several models. 
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Figure 1: a-Isobars and lines of maximum attainable Ra in the density-temper at urc plane, 
b- Semi log enlargement of the plot a, corresponding to the experiments described in the 
text. The maximum Ra achieved within Boussinesq conditions in the 2 cm high cell is 10^^. 
The dark area corresponds to the 2-phases region. 

Lines of maximum attainable Rayleigh numbers, hereafter called iso-Rttmax curves, 
have been calculated at various pressures and temperatures for the mini-cell and for the 
maximum temperature difference ATj^^x compatible with a Boussinesq criterion [2], here 
defined as: a AT < 20%. This criterion is close to the one respected in the experimental 
data for which a AT < 21%. In density-temperature coordinates iso-Ramax curves arc 
shown on figures 1-a and 1-b as the thick full lines. Close to the critical point, the divergence 
of ^, easily accessible in ^He (2.2 bars, 5.2 K) is clearly illustrated by the extremely high 
Ra which can be obtained in reasonable experimental conditions (AT > 1 mK). The thin 
lines on figure 1-a correspond to isobars. 

In the large cells[6, 11], the Pr variation is only obtained for Ra above 10^°. Clear 
understanding of the Pr variation is difficult in these cells due to the occurrence of the ul- 
timate regime. With cell dimensions divided by ten, the Pr variation is already observable 
around Ra — lO"^. 

It is worth noticing that the divergence in Ra^ax is due to Cp in the relation ^ = ^x~' 
77 being the fluid viscosity and Cp the speciflc heat at constant pressure. The Prandtl 
number varies as: Pr = . The Cp divergence appears far away from the critical point 
and gives a long range effect to the rapid variation of Ra and Pr. In various experimental 
conditions, AT has been widely varied up to three decades at roughly constant temperature 
and density. This is the largest excursion ever achieved in Rayleigh-Benard experiments 
within Boussinesq condition. This gives access to power law exponent Nu versus Ra at 
constant Pr and independently of the fluid properties knowledge. Indeed, as temperature 
and density in the bulk of the flow are almost constant, the same hold for the He properties^. 
Thus, the power law exponent of Nu{Ra) is independent of possible error on the fluid 
properties, if we except the adiabatic gradient correction. However these properties have 
to be known precisely for the Nu versus Pr dependence studies. 

"'^In our measurements, the average temperature and density slightly differ from one point to another. 
These variations in the experimental conditions are precisely measured and the fluid properties are recalcu- 
lated for each point. If we make the unrealistic hypothesis that the fluid properties variation is estimated 
with a 100% error, the resulting uncertainty of the effective power law of the Nu(Ra) dependence would 
be less than 4%. 
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Figure 2: Sclieme of tlie experimental set-up. 



When this experiment was designed, the situation was the following : two experiments [12, 
13] were given contradictory results, each being is agreement with a different theory. The 
first one, conducted over nearly 2 decades of Pr (1 < Pr < 93) found a —0.2 exponent 
for the effective Nu{Pr) power law, while the second experiment's data can be fitted with 
a -0.01 exponent over 0.9 decade of Pr (4 < Pr < 34). Our aim was to elucidate the 
controversy and to expand the explored Pr range below Pr = 1. Our experiment allows 
to vary the Prandtl over 1.5 decade, that is 0.7 < Pr < 21 and for a large excursion of Ra 
numbers. The references [14, 15, 16, 17] present Pr-dependence studies conducted in the 
Rayleigh-Benard geometry for lower Ra, much lower or higher Pr, or for a much smaller 
Pr range. The thermal control and measurement accuracy of our experiment are unprece- 
dented in cryogenics convection experiments for Ra > 3.10^. It revealed two unexpected 
effects : the side-wall effect and the bimodality. 

This paper gives a detailed description of the apparatus (section 2) and improvements 
of ^He properties fits (section 3). In section 4, we remind the three main results : the 
bimodality of the flow, the side-wall effect and the Prandtl number dependence. Section 5 
proposes some perspectives for convection studies. 



2 INSTRUMENTATION 

2.1 The Rayleigh-Benard cell design 

The experimental set-up, presented on figure 2, is placed in a cryogenic vacuum. The 1/2 
aspect ratio mini-cell {h = 2 cm), is also shown on figure 3. The stainless-steel cylindrical 
wall is 0.25 mm thick. It can hold pressures up to several tens of bars. The upper plate is 
part of a main Cu flange which ensures the thermal link to the liquid He bath, through a 
brass plate of measured thermal resistance (53 K/W at 4.3K) and a high conductivity Cu 
post [18, 19, 20, 21]. The brass plate acts as a thermal resistance which allows to regulate 
the top plate at temperatures different from that of the He bath. The lower plate is also 
made of copper. In such a helium/copper set-up and for the Ra numbers explored in 
this study, the plates properties (finite conductivity and heat capacity) do not alter the 
dynamical formation of the coherent structures (plumes,...) [22, 23]. 

Special care has been taken during the cell assembly. The roughness of the plates 
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Figure 3: 4 a- The Rayleigh-Benard cell. 4b- Photograph of the experimental set up (same 
arrangement as on fig. 2 ). 

surfaces in contact with the fiuid is estimated less than 2 fim. The silver soldering of the 
wall and the copper rings ensure the connection to the plates. For all the cells developed 
in our group, the side wall design is chosen in order to have a perfect cylinder along all 
the active length of the cell[18, 19, 20, 21]. The plates parallelism is guaranteed by special 
machining procedure performed after the silver soldering: the distance between plates is 
19.99 ±0.02 mm. 

The cell is filled through a capillary closed with a cold needle valve located in the main 
helium bath. This capillary is connected in series with a capacitance cell, that will be 
further described, and the Rayleigh-Benard cell. 



2.2 Experimental procedure 

Each cell plate holds two Ge resistors from the same batch with close resistance values 
and temperature dependence. The upper fiange is temperature regulated with a PID 
(Proportional-Integral-Derivative) analogue regulator with a fifth Ge resistor. The tem- 
perature difference AT between the plates in the Rayleigh-Benard cell is determined from 
the measurement of the ratio of two Ge resistors using a resistance ratio bridge operating 
at 30 Hz with 1 /lA current amplitude [24]. The resistance ratio variation with and without 
heating gives AT through the calibration of the resistors and the additional measurement 
of the upper plate temperature. This procedure is valid even for AT larger than 1 K, as 
checked with the direct temperature measurement of each plate. The ratio without heating 
(zero AT) is monitored during twelve hours before and after each measurement cycle. In 
less than one hour the equilibrium value is obtained, except for the data close to the critical 
point: in such conditions the thermal diffusion time diverges and the used zero AT is the 
one obtained at a lower density. 

A ratio variation of 10~^ at 5 corresponds to 25 fiK for AT. The stability of the 
set-up and electronic apparatus is better than 30 fiK over 12 hours. The radiation heat 
losses are estimated to be around 10 nW which gives typically 50 fiK for AT. This is 
of the order of the adiabatic gradient temperature difference [2, 20] in the 2 cm high cell, 
this effect limiting the smallest achievable AT. More than 3 decades of variation of AT 
have been achieved for given mean temperature and density, from below the convection 
onset up to the turbulent regimes. The Boussinesq criterion is a. AT < 21%. Besides the 
conductivity and viscosity never vary by more than 8% between the bottom and top of 
the cell. Still between the bottom and top of the cell, the constant-pressure heat capacity 
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and the expansion coefficient a vary by less tfian 10% for 88% of tfie points, they vary 
by less than 20% for 95% of the points and by less than 40% for all the points. In large 
cells, an original thermocouple technique is more appropriate than a resistance bridge to 
measure the temperature difference: the AT zeroing procedure is not compatible with the 
large thermal relaxation times. 

A four wires voltage standard (Electronics Development Corporation) provides a con- 
stant heating power on the lower plate. In order to limit local overheating, the heater 
is distributed on the surface. All the copper pieces are made out of a commercial non- 
annealed Cu, which was characterized in another experiment: its thermal conductivity is 
around 400 W jm.K. We have measured the wall thermal conductance between 4.5 and 
6 K. For cross-validation, two measurements have been done with an empty cell and with 
helium at 80 gjw?. After subtraction of the diffusive helium contribution both results 
agree within 2.5%: in this difference 2% are explained by the cell design[21]. This side-wall 
contribution is described by the fit —40.1 + 44.75.T in [iW jK including the conductance 
of the copper heating wires (33 puWj K at 4.5 K with 20% uncertainty) going to the lower 
plate. The lower plate heat capacity, as measured by a relaxation method, is 73 mJ/K 
including addenda (Ge resistor holders and copper ring). The brass plate heat leak has a 
measured resistance of 53 K/W at 4.3 K. 



2.3 Densitometry 

In order to determine the Ra and Nu with an absolute resolution of a few percent, a density 
accuracy of 1%, at least, is needed. The dead volume coming from the filling capillary 
going to room temperature is too large to determine with enough precision the He density 
during the cell filling procedure. We thus have performed an in-situ measurement using a 
capacitive probe located in a specific cell, in order to have a much better resolution. The 

density is extracted from the Clausius-Mosotti relation: ^ = where e, p, M, Q are 
respectively the permittivity, the density, the molar mass and the polarisability of helium 
(C = 0.123296 cmVmoZ) [25]. 

Two capacitances are placed in a ratio bridge. The two porous frames of the "ac- 
tive" capacitance (C ~ 17.5 pF) are made out of printed circuits and 0.1mm separated. 
This capacitance is located in the capacitive cell and totally immersed in helium. Under 
these conditions the mechanical dependence with pressure effects is minimized. On the 
inner part of each frame a circular electrode (16 mm diameter) is engraved together with 
a guard ring. Special attention in the design reduces differential contraction effect and 
parasitic capacitances : no spurious effect were detected and no temperature effects were 
observed. The other capacitance (C^ = 9.4 pF) made out of mica, is located at 4.2 K in 
the calorimeter vacuum. It is the reference one. 

The bridge operates at 3 kHz. The ratio between both capacitances is a direct mea- 
surement of e. The density measurement range is — 140 kg/m^, under pressures from 
up to 7 bars and temperatures between 4.5 and 6.5 K. Over ten days the stability is 
10"^ (±40gf/m^), that is better than 0.1% in density for Ra above 2.10^. Two calibrations 
of the capacitance ratio at the beginning and at the end of the experiment agree within 
30 g/m^. The signal averaged over 30 s has a resolution of 10~^ (less than 1 g/m^), which 
can be maintained over a few hours. 

In principle k = should be a constant over the density range. We have achieved low 

and high densities calibrations. During the low pressure calibration, the cell, connected to a 
few litres reservoir at room temperature, is regulated at 5.432 K. For pressures lower than 
1000 mbar, no condensation occurs in the filling line. In order to have stable operation 
conditions we restrain the low-pressure calibration below 350 mbars: a precise pressure 
measurement gives access to the density through ref. [26] . At high density we measured 
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Figure 4: Calibration of the density measurement cell. Insert: Difference, in percent, 
between the density estimations with and without pressure correction. 

the density following a procedure described below for the absolute temperature calibration 
(see Thermometry). For example, with 31 /iW heating power applied on the bottom plate, 
we extract a density of 112.73 kg/m? from the measured boiling temperature (4.7088 K 
with AT = Q.hmK). 

The calibration results are summarized on figure 4 where is plotted as a function 
of the pressure P. varies by 6.10"^ over the whole range and we have assumed a linear 

variation versus P. Such behaviour is typical of a residual mechanical deformation. On 
the figure insert the difference between density with and without the linear correction is 
plotted versus p. This 5p is less than 0.5% and goes through a minimum at 0.15% close 
to the critical density. In all the following we obtain the density from the linear pressure 
correction and we estimate the p uncertainty to be about 0.1%. 

When both cells are filled and the needle valve is closed, no temperature dependence of 
the capacitive signal is expected. However a 0.3% tiny reproducible (on a few months scale) 
variation was observed as seen on figure 5. This was explained by the helium compression in 
the upper part of the filling line close to the needle valve in thermal contact with the main 
helium bath at 4.22 K. We have evidenced a linear correlation through the comparison of 
the total measured density and the calculated density in the capillary. The slope, plotted 
on the insert of figure 5 is the ratio between the total volume and the capillary dead 
volume : the value 23 ±4 is in fair agreement with a less precise value extracted from 
a geometric determination. The ±4 uncertainty on the volume ratio (due to the scatter 
of points) corresponds to a measured density uncertainty less than 0.06% on the whole 
densities and temperature range. It confirms by an independent way our formerly quoted 
0.1% density uncertainty. Note that this capillary effect has strictly no influence on the 
density measurement in the cells. 



2.4 Thermometry 

A one millikelvin uncertainty on the mean temperature gives an uncertainty up to 1% on 
the Ra and Nu values, in the range of temperature and density of this experiment. In order 
to compare the various data issued from several references [25, 26, 27, 28, 29, 30, 31, 32, 33], 
we use the ITS-90[34] critical temperature: = 5.1954 K and adjust the thermometer 
in-situ calibration onto that value. This calibration procedure is illustrated on figure 6 at 
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Figure 5: Measurement of the helium compression effect in the filhng capillary. Insert : 
Measured ratio of the cells volume divided by the 4.22 K capillary volume versus density. 
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Figure 6: Full squares: Cell density versus top plate temperature around the condensation 
onset. Dotted circles: Plates temperature difference in arbitrary units. 



a density of 51.5 kg/m? where the condensation is expected at: T = Tc — 47.4 mK. We 
apply a small 125 fiW heating power and monitor both AT across the Rayleigh-Benard 
cell and the gas density pgas in the capacitive cell. The temperature of the upper plate 
is slowly lowered. The sharp drop of AT and pgas is the signature of the condensation in 
the cell. It is worth to note that the density measurement is more precise than the AT 
one [12, 19] to identify condensation and enables a 1 mK resolution. 



3 ^He PROPERTIES 



3.1 Thermal expansion coefficient 



The reference fits of Arp and McCarty[30, 26] account for ^He thermodynamics properties 
over a wide temperature and pressure range (0.8 — 15000 K, — 2000 MPa) but ignores 
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Figure 7: Comparison between the thermal expansion coefficients at density p = 
{±0.6kg/m^ for the experimental data). Present measurements are in full circles. 
Kierstead[29] (T < 5AK) and McCarty and Arp[26] (T > 5AK) fits (lines) have been 
extrapolated in the circled area. The XHePak fit (pluses) is the commercial fitting package 
cited in ref. [7] 



the critical divergence, which has been fitted by Kierstead[28, 29]. Unfortunately, the tem- 
perature and density validity ranges of the two fits do not overlap. Besides, extrapolation 
of the thermal expansion coefficient from both ffis suggests up to 25% discrepancy between 
these ffis, which is incompatible with the extrapolation uncertainty. Figure 7 illustrates 
this discrepancy for a density of p = lAkg/m?. 

In our experiment, the high sensitivity of the capacitive cell gives access to a, in a 
temperature and pressure range overlapping both ffis. The Rayleigh-Benard cell heating 
increases the density in the capacitive cell. For low AT across the convection cell, the 
density variation is given by: 5p = ^pa ^^^^^ , where vi and V2 are the capacitive cell and 

Rayleigh-Benard cell volumes respectively. The results are shown on figure 8. We extract 
the slope of these curves for AT going to 0, for given temperature and density conditions. 
We have also done a correctional] to take into account the volume of the 4.2 K ffiling 
capillary below the needle valve: this correction represents a few percent at low and high 
densities and is smaller than 0.5% between 40 and 80 kg/m?. Each point corresponds to 
a Temperature-Density condition which fully falls into the validity range of either one of 
the two ffis. 

Determination of the geometrical coefficient ^^^^^ is illustrated on figure 9 (the capillary 

volume introduces less than 1% correction on this formula and this correction is not shown 
here but it is taken into account in the analysis of ref. [21]). Figure 9 shows the quantity 

Qest/ ^"^!^^ = 01 est I wherc aest is an estimated a value from the literature [26, 29, 30]. 
The X-axis is chosen in order to avoid the degeneracy for the data taken at same density 
but at various temperatures. Most of the data are compatible with the value 11.15 ±0.25. 

As expected the data using the values from ref. [26] are very reliable far from the 
critical point, but they need a correction which rises up to about 20% when approaching 
this regime. The data from ref. [30] are in good agreement with those of ref. [26] at low 
density but need a correction of several percent at high density. In the critical region 
ref. [29] appears to give the best agreement, as illustrated on Figure 7. The various a 
values are summarized in Table 1. 
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Figure 8: Capacitive measurement of the p increase in the capacitive cell due to the 
Rayleigh-Benard cell heating versus the plates temperature difference in the convection 
cell. 
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Figure 9: Geometrical coefficient (see text) versus p + 5T for the determination of the 
helium expansion coefficient. 
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Table 1 : 



ocmeasur. '■ expansion coefficient measurements, Aameasur. total uncertainty on ameasur., 
OiMCA '■ expansion coefficient estimated with McCarty and Arp 1990 fit [26] ; OiKierstead '■ 
expansion coefficient estimated with Kierstead fit [29]. 

In our Rayleigh-Bcnard data analysis, the a coefficient and other thermodynamical 
coefficients such as {Cp — C„) are obtained from the first order partial derivatives of the 
pressure versus temperature or density. Such a way of deriving properties ensures the self- 
consistency between thermodynamics parameters. As a consequence, taking into account 
our measurements we have corrected the fits at their source, that is directly on the fit of 
{dP/dp)T,MCA ofref. [26] : 

(dP/dp)T = {dP/dp)T,MCA[l + Fr.Fp] 

with Ft = 4.62 - 0.6587 and Fp = 0.246 - 0.001 17(p - 67)^ 

Thus our recommendations are the following : in the critical zone as defined by Kier- 
stead (55.7 < p < 83.5 kg/m^' and T < 5.362 K) use the Kierstead fit, out of this range, 
use the Arp and MacCarty[26] fit with the above correction if both Ft and Fp are pos- 
itive and if the total relative correction FT-Fp is larger than 3%. With this correction, 
Kierstead fit (unchanged) and Arp and McCarty fit (modified and extrapolated) reconnect 
much better. 

We should mention here that Cy cannot be derived from the state equation. In the 
zero-density region, is known exactly (perfect gas) and in the critical region, it has 
been fitted by Moldover[27]. In between, we resorted to exact thermodynamics relation to 
bridge to either one of these two regions. 



3.2 The transport properties : viscosity and thermal conductivity 

A fit of "^He viscosity in the range 4 — 20 and — 10 MPa has been proposed by Stew- 
ard and Wallace[35]. In our range of parameters, the fit is an interpolation of isothermal 
measurements at 4, 5, 6 and 10 K conducted by these authors. Along the critical iso- 
chore, comparison with viscosity data of Kogan et a/. [36] and Agosta et a/. [33] shows +7% 
deviation at 5.2 K and —7% at 7 K. However Steward and Wallace measurements at 4, 
5 and 6 K are in a few percent agreement with the literature, including the 2 references 
mentioned above, but their data at 10 K differ significantly from the literature. It ap- 
peared that this 10 isothermal entails a strong bias on Steward and Wallace fit down 
to the lower temperatures : this is consistent with a concern regarding a contamination 
of helium, due to a defective purifier[37]. Consequently, we derived a new interpolation 
between the 4, 5 and 6 K isothermals above 70 kg/m^ and with additional data along the 
critical isochore[36], in the zero-density limit (ab-initio calculation of ref. [38]), and on the 
vapour-liquid curve[25]. Concerning viscosity, there is a clear need of new measurements 
in the range 6K — lOK and above p^. We have no data in this range but this lack of 
information makes doubtful the interest in publishing our fit. 

The thermal conductivy has been estimated from a specially designed new fit through the 
data of Acton and Kellner[31, 32]. We re-computed the density data of these papers, which 
had been estimated with the 1973 fit of McCarty [30], even in the critical region. Our new 
fit agrees within ±2% with the published [31, 32] and unpublished [36] data of Acton and 
Kellner. Whenever it was possible, our convection measurements have been conducted 
at the same mean temperatures as the one employed by Acton and Kellner, in order to 
minimize interpolation errors. 



11 



4 EXPERIMENTAL RESULTS AND ANALYSIS 



The remaining of this paper presents the heat transfer measurements and their conse- 
quences. First we consider the Nu{Ra) relation. Our high accuracy on Nu gives access to 
two new effects : the bimodahty of Nu and the side-wall conductivity influence on heat 
transfer. We then turn to the Nu{Pr) dependence. 



4.1 Bimodality 

On figure 10 we show Nu as a function of Ra. In order to display all data with an increased 
vertical resolution, Nu is arbitrarily re-normalized by i?a°'^^. For comparison, we display 
the data at Pr = 0.7 and 1.1 from ref. [8]. The change in slope at Ra ~ 10^ — 2.10^ 
corresponds to the soft-hard turbulent transition in 1/2 aspect ratio cell[39]. In the soft 
regime, the exponent of Nu vs. Ra is close to 0.25 (interpolation over only 1.5 decade), 
while in the hard turbulence regime the exponent is close to 0.31, in between the 2/7 and 1/3 
predictions of traditional theories [40]. It is interesting to note that the exponent averaged 
over these two regimes is close to 2/7. Another possible interpretation of the exponent 
change would be to reject the soft-hard transition picture and rather see a continuous 
variation of the exponent resulting from a linear combinaison of two power laws, although 
the abruptness of the exponent change doesn't fully fit in this picture. 

For 2.10^ < Ra < 2.10^°, points can be gathered into two subsets which differ by roughly 
5 — 7% in Nu. Such a data bimodality cannot be taken into account by the uncertainties, 
which arc twice smaller than this Nu gap, nor by a Pr dependence. Switches from one set of 
data to the other occur varying AT under quasi- const ant mean temperature and density 
conditions (see for example Pr ~ 0.95) : this definitely rules out that the bimodality 
would come from an improper helium property estimation in the T-p plane. A numerical 
simulation conducted for the same cell geometry[41] recently found that two types of large 
scale fiow can fit the cell. This mechanism of bimodality is consistent with the invariance 
of the Nu gap (in log scale) observed in our data. 

We cannot decide if the bimodahty reveals spurious effect of the boundary conditions or 

a macroscopic degree of freedom of the fiow with a slow dynamics. In the first hypothesis, 
each mode could be stabilized by the thermal inertia of cell boundaries (for instance, 
an ascending wind warms up the nearby side-wall which -in turns- enhance an ascending 
convection) and the switching from one mode to the other should be hysteretic. We arc not 
able to precise more what is the anchoring mechanism. The second hypothesis has drastic 
consequences since the very slow dynamics (at least hundreds of turn-over times) ruins -for 
practical reasons- the present definition of Nu : indeed a clean averaging procedure would 
request a time duration incompatible with a laboratory experiment. 

In this paragraph, we present a practical difference between the cryogenic and room 
temperature experiments. These latter ones are illustrated quantitatively with water as 
the working fluid. The characteristic time of convection is : 



At a given Ra and Pr, the caracteristic time tne and t^jater in helium and water obeys 
to : 




(6) 
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p 


T 


P 


Q 


AT 


ATadiab 


a.AT 


w 


Pr 


Ra 


Nu 




w 




















6.963 


4.489 


0.566 


3.326 


4.681 


44 


0.1 4 


3.65 


0.84 


3.049 E 6 


1 1 .4 


6.961 


4.490 


0.566 


5.102 


6.557 


44 


0.20 


3.65 


0.84 


4.277 E 6 


1 2.8 


6.954 


4.491 


0.566 


7.801 


9.353 


44 


0.29 


3.65 


0.84 


6.094 E 6 


1 3.9 


6.954 


4.491 


0.566 


7.801 


9.353 


44 


0.29 


3.65 


0.84 


b.094 E b 


1 3.9 


6.955 


4.493 


0.566 


1 2.05 


1 3.1 9 


44 


0.41 


3.65 


0.84 


8.591 E 6 


1 5.5 


6.947 


4.496 


0.566 


1 8.56 


1 8.67 


44 


0.57 


3.65 


0.84 


1 .21 2 E 7 


1 7.1 


6.947 


4.496 


0.566 


1 8.56 


1 8.67 


44 


0.57 


3.65 


0.84 


1 .21 2 E 7 


1 7.1 


6.779 


4.557 


0.563 


28.02 


26.38 


44 


0.79 


3.67 


0.83 


1 .542 E 7 


1 8.2 


6.779 


4.557 


0.563 


28.02 


26.38 


44 


0.79 


3.67 


0.83 


1 .542 E 7 


1 8.2 


6.779 


4.557 


0.563 


28.02 


26.58 


44 


0.79 


3.67 


0.83 


1 .554 E 7 


1 8.1 


6.779 


4.557 


0.563 


28.02 


26.58 


44 


0.79 


3.67 


0.83 


1 .554 E 7 


1 8.1 


6.944 


4.499 


0.566 


25.59 


24.66 


44 


0.76 


3.65 


0.84 


1 .595 E 7 


1 8.0 


6.774 


4.556 


0.563 


34.47 


30.98 


44 


0.93 


3.67 


0.83 


1 .809 E 7 


1 9.3 


6.774 


4.556 


0.563 


34.47 


30.98 


44 


0.93 


3.67 


0.83 


1 .809 E 7 


1 9.3 


6.773 


4.557 


0.563 


34.44 


31 .88 


44 


0.95 


3.67 


0.83 


1 .861 E 7 


1 8.6 


6.773 


4.557 


0.563 


34.44 


31.88 


44 


0.95 


3.67 


0.83 


1 .861 E 7 


1 8.6 


6.933 


4.502 


0.566 


35.49 


32.03 


44 


0.98 


3.66 


0.84 


2.060 E 7 


1 9.5 


6.933 


4.502 


0.566 


35.49 


32.03 


44 


0.98 


3.66 


0.84 


2.060 E 7 


1 9.5 


6.768 


4.556 


0.562 


42.34 


37.37 


44 


1 .1 


3.67 


0.83 


2.1 78 E 7 


1 9.7 


6.768 


4.556 


0.562 


42.34 


37.37 


44 


1 .1 


3.67 


0.83 


2.1 78 E 7 


1 9.7 


6.767 


4.556 


0.562 


42.33 


37.67 


44 


1 .1 


3.67 


0.83 


2.1 96 E 7 


1 9.5 


6.767 


4.556 


0.562 


42.33 


37.67 


44 


1 .1 


3.67 


0.83 


2.1 96 E 7 


1 9.5 


6.761 


4.556 


0.562 


52.59 


43.87 


44 


1 .3 


3.68 


0.83 


2.552 E 7 


21 .1 


6.761 


4.556 


0.562 


52.59 


43.87 


44 


1 .3 


3.68 


0.83 


2.552 E 7 


21.1 


6.761 


4.556 


0.562 


52.58 


44.07 


44 


1 .3 


3.67 


0.83 


2.564 E 7 


21 .0 


6.761 


4.556 


0.562 


52.58 


44.07 


44 


1 .3 


3.67 


0.83 


2.564 E 7 


21 .0 


6.927 


4.507 


0.567 


49.03 


42.06 


44 


1 .3 


3.66 


0.84 


2.690 E 7 


20.7 


6.754 


4.555 


0.561 


64.95 


50.96 


44 


1 .5 


3.68 


0.83 


2.959 E 7 


22.7 


6.754 


4.555 


0.561 


64.95 


50.96 


44 


1 .5 


3.68 


0.83 


2.959 E 7 


22.7 


6.752 


4.556 


0.561 


64.90 


52.56 


44 


1 .6 


3.68 


0.83 


3.049 E 7 


21 .9 


6.752 


4.556 


0.561 


64.90 


52.56 


44 


1 .6 


3.68 


0.83 


3.049 E 7 


21 .9 


6.91 


4.51 4 


0.566 


67.79 


54.84 


44 


1 .7 


3.66 


0.84 


3.471 E 7 


22.1 


6.91 


4.51 4 


0.566 


67.79 


54.84 


44 


1 .7 


3.66 


0.84 


3.471 E 7 


22.1 


6.742 


4.556 


0.560 


80.12 


62.65 


44 


1.9 


3.68 


0.83 


3.620 E 7 


22.8 


6.742 


4.556 


0.560 


80.1 2 


62.65 


44 


1 .9 


3.68 


0.83 


3.620 E 7 


22.8 


6.741 


4.557 


0.560 


80.1 2 


63.05 


44 


1 .9 


3.68 


0.83 


3.642 E 7 


22.6 


6.741 


4.557 


0.560 


80.1 2 


63.05 


44 


1 .9 


3.68 


0.83 


3.642 E 7 


22.6 


6.906 


4.51 8 


0.567 


81.43 


63.96 


44 


1 .9 


3.66 


0.84 


4.029 E 7 


22.9 


6.731 


4.557 


0.560 


98.04 


73.04 


44 


2.2 


3.68 


0.83 


4.204 E 7 


24.2 


6.731 


4.557 


0.560 


98.04 


73.04 


44 


2.2 


3.68 


0.83 


4.204 E 7 


24.2 


6.729 


4.558 


0.560 


97.98 


75.04 


44 


2.2 


3.68 


0.83 


4.31 3 E 7 


23.4 


6.729 


4.558 


0.560 


97.98 


75.04 


44 


2.2 


3.68 


0.83 


4.31 3 E 7 


23.4 


6.893 


4.523 


0.566 


98.01 


73.90 


44 


2.2 


3.66 


0.84 


4.61 6 E 7 


24.0 


6.893 


4.523 


0.566 


98.01 


73.90 


44 


2.2 


3.66 


0.84 


4.61 6 E 7 


24.0 


6.71 6 


4.559 


0.559 


1 21 .4 


88.33 


43 


2.6 


3.68 


0.83 


5.050 E 7 


24.9 


6.71 6 


4.559 


0.559 


1 21 .4 


88.33 


43 


2.6 


3.68 


0.83 


5.050 E 7 


24.9 


6.71 6 


4.558 


0.559 


1 21 .4 


88.33 


43 


2.6 


3.68 


0.83 


5.051 E 7 


24.9 


6.71 6 


4.558 


0.559 


1 21 .4 


88.33 


43 


2.6 


3.68 


0.83 


5.051 E 7 


24.9 


6.886 


4.529 


0.567 


1 17.7 


85.94 


44 


2.6 


3.66 


0.84 


5.332 E 7 


24.9 


6.699 


4.560 


0.558 


1 51 .2 


1 05.2 


43 


3.1 


3.68 


0.83 


5.972 E 7 


26.2 


6.699 


4.560 


0.558 


1 51 .2 


1 05.2 


43 


3.1 


3.68 


0.83 


5.972 E 7 


26.2 


6.698 


4.561 


0.558 


1 51 .1 


105.6 


43 


3.1 


3.68 


0.83 


5.992 E 7 


26.1 


6.698 


4.561 


0.558 


1 51 .1 


1 05.6 


43 


3.1 


3.68 


0.83 


5.992 E 7 


26.1 


6.872 


4.536 


0.567 


1 41 .6 


1 00.1 


44 


3.0 


3.67 


0.84 


6.1 50 E 7 


25.9 


6.680 


4.563 


0.557 


184.0 


124.1 


43 


3.7 


3.68 


0.83 


6.982 E 7 


27.2 


6.680 


4.563 


0.557 


184.0 


124.1 


43 


3.7 


3.68 


0.83 


6.982 E 7 


27.2 


6.679 


4.564 


0.557 


183.9 


124.3 


43 


3.7 


3.68 


0.83 


6.989 E 7 


27.1 


6.679 


4.564 


0.557 


183.9 


124.3 


43 


3.7 


3.68 


0.83 


6.989 E 7 


27.1 


6.670 


4.565 


0.557 


201.6 


133.8 


43 


4.0 


3.68 


0.83 


7.490 E 7 


27.7 


6.670 


4.565 


0.557 


201.6 


133.8 


43 


4.0 


3.68 


0.83 


7.490 E 7 


27.7 


6.669 


4.566 


0.557 


201.5 


135.1 


43 


4.0 


3.68 


0.83 


7.556 E 7 


27.4 


6.669 


4.566 


0.557 


201.5 


135.1 


43 


4.0 


3.68 


0.83 


7.556 E 7 


27.4 



with obvious notations. In helium and water cells, the smallest convection times are 
obtained at the highest a.AT which -given Boussinesq conditions- are typically (a. AT)^e — 
0.2 and (a. AT)^,^^^ — 0.02. The Ra explored in this hue = 2 cm experiment are achievable 
with water for h^jater = 20 cm typically, which gives twater/tne = 10. We can state more 
generally that the times scales in cryogenic helium are typically 10 times smaller than in 
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p 


T 


P 


Q 


AT 


ATadiab 


a.AT 


w 


Pr 


Ra 


Nu 


IKg/m-'j 
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6.656 


4.567 


0.556 


228.0 


1 47.5 


43 


4.4 


3.68 


0.83 


0.2U2 E 7 


28.6 


6.656 


4.567 


0.556 


228.0 


1 47.5 


43 


4.4 


3.68 


0.83 


O.C.QZ E 7 


28.6 


6.654 


4.569 


0.556 


227.9 


1 49.9 


43 


4.4 


3.68 


0.83 


8.321 E 7 


28.0 


6.654 


4.569 


0.556 


227.9 


149.9 


43 


4.4 


3.68 


0.83 


8.321 E 7 


28.0 


6.643 


4.570 


0.555 


253.0 


1 60.5 


43 


4.7 


3.68 


0.83 


8.873 E 7 


29.2 


6.643 


4.570 


0.555 


253.0 


1 60.5 


43 


4.7 


3.68 


0.83 


8.873 E 7 


29.2 


6.642 


4.570 


0.555 


253.0 


1 61 .3 


43 


4.8 


3.68 


0.83 


8.91 2 E 7 


29.0 


6.642 


4.570 


0.555 


253.0 


1 61 .3 


43 


4.8 


3.68 


0.83 


8.91 2 E 7 


29.0 


6.630 


4.573 


0.555 


280.6 


1 73.8 


43 


5.1 


3.68 


0.83 


9.541 E 7 


30.0 


6.630 


4.573 


0.555 


280.6 


1 73.8 


43 


5.1 


3.68 


0.83 


9.541 E 7 


30.0 


6.628 


4.573 


0.555 


280.5 


1 75.7 


43 


5.2 


3.68 


0.83 


9.636 E / 


29.6 


6.628 


4.573 


0.555 


280.5 


1 75.7 


43 


5.2 


3.68 


0.83 


9.636 E 7 


29.6 


6.61 3 


4.576 


0.554 


31 1 .5 


1 90.1 


43 


5.6 


3.68 


0.83 


1 .035 E 8 


30.5 


6.61 3 


4.576 


0.554 


31 1 .5 


1 90.1 


43 


5.6 


3.68 


0.83 


1 .035 E 8 


30.5 


6.61 3 


4.576 


0.554 


31 1 .5 


1 91 .0 


43 


5.6 


3.68 


0.83 


1 .039 E 8 


30.4 


6.61 3 


4.576 


0.554 


31 1 .5 


1 91 .0 


43 


5.6 


3.68 


0.83 


1 .039 E 8 


30.4 


6.610 


4.577 


0.554 


319.0 


194.0 


43 


5.7 


3.68 


0.83 


1 .054 E 8 


30.6 


6.61 


4.577 


0.554 


319.0 


1 94.0 


43 


5.7 


3.68 


0.83 


1 .054 E 8 


30.6 


6.608 


4.578 


0.554 


318.9 


1 95.4 


43 


5.8 


3.68 


0.83 


1 .060 E 8 


30.4 


6.608 


4.578 


0.554 


318.9 


1 95.4 


43 


5.8 


3.68 


0.83 


1 .060 E 8 


30.4 


6.573 


4.585 


0.552 


395.6 


231 .0 


43 


6.8 


3.69 


0.83 


1 .231 E 8 


32.1 


6.573 


4.585 


0.552 


395.6 


231.0 


43 


6.8 


3.69 


0.83 


1 .231 E 8 


32.1 


6.572 


4.585 


0.552 


395.5 


231 .9 


43 


6.8 


3.69 


0.83 


1 .235 E 8 


32.0 


6.572 


4.585 


0.552 


395.5 


231.9 


43 


6.8 


3.69 


0.83 


1.235 E 8 


32.0 


6.529 


4.596 


0.551 


491.7 


275.4 


43 


8.0 


3.69 


0.82 


1 .434 E 8 


33.6 


6.529 


4.596 


0.551 


491 .7 


275.4 


43 


8.0 


3.69 


0.82 


1 .434 E 8 


33.6 


6.528 


4.596 


0.551 


491 .6 


276.6 


43 


8.1 


3.69 


0.82 


1 .439 E 8 


33.5 


6.528 


4.596 


0.551 


491 .6 


276.6 


43 


8.1 


3.69 


0.82 


1 .439 E 8 


33.5 


6.477 


4.61 


0.549 


61 1 .7 


329.5 


43 


9.5 


3.70 


0.82 


1 .665 E 8 


35.1 


6.477 


4.61 


0.549 


61 1 .7 


329.5 


43 


9.5 


3,70 


0.82 


1 .665 E 8 


35.1 


6.476 


4.61 


0.549 


61 1 .7 


330.1 


43 


9.6 


3.70 


0.82 


1 .668 E 8 


35.0 


6.476 


4.61 


0.549 


61 1 .7 


330.1 


43 


9.6 


3.70 


0.82 


1 .668 E 8 


35.0 


9.699 


5.500 


0.959 


1 3.52 


1 1 .82 


45 


0.30 


3.87 


0.84 


8.6Z7 E 


1 5.7 


9.693 


5.501 


0.959 


1 9.51 


16.05 


45 


0.41 


3.88 


0.84 


1 .1 71 E 7 


16.9 


9.687 


5.500 


0.958 


28.1 7 


21 .66 


45 


0.56 


3.88 


0.84 


1 .579 E 7 


1 8.3 


9.680 


5.501 


0.958 


40.63 


29.53 


45 


0.76 


3.88 


0.84 


2.149 E 7 


1 9.6 


9.669 


5.499 


0.957 


58.61 


39.78 


45 


1.0 


3.88 


0.84 


2.892 E 7 


21.3 


9.655 


5.500 


0.956 


84.57 


54.14 


45 


1.4 


3.88 


0.83 


3.921 E 7 


22.9 


9.636 


5.499 


0.954 


1 21 .2 


72.75 


45 


1.9 


3.88 


0.83 


5.247 E 7 


24.7 


9.608 


5.498 


0.951 


1 74.8 


98.07 


45 


2.5 


3.88 


0.83 


7.029 E 7 


26.8 


9.573 


5.499 


0.949 


253.9 


1 31 .3 


44 


3.4 


3.88 


0.83 


9.325 E 7 


29.5 


9.525 


5.500 


0.945 


365.8 


1 79.0 


44 


4.6 


3.88 


0.83 


1 .255 E 8 


31 .6 


9.461 


5.499 


0.939 


527.4 


240.9 


44 


6.2 


3.88 


0.83 


1 .664 E 8 


34.3 


9.372 


5.502 


0.932 


760.4 


324.2 
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8.2 


3.88 


0.83 


i2.1 8b E 8 


37.2 


9.260 


5.500 


0.922 
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434.1 
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1 1 


3.89 


0.83 


il.o4D E O 


40.6 


9.859 


4.808 


0.824 


2.1 53 


2.584 


46 
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3.68 


0.90 


3.078 E 6 


1 2.7 


9.856 


4.809 


0.824 


4.475 


4.787 


46 
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3.68 


0.90 


5.741 E 6 


1 4.5 


9.854 


4.809 


0.824 


6.488 


6.560 


46 


0.21 


3.68 


0.90 


"7 OOf^ C f 

l.oob E 6 


1 5.5 


9.851 


4.806 


0.823 


9.390 


8.929 


46 


0.28 


3.68 
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1 .077 E 7 


1 6.6 


9.847 


4.806 


0.823 


1 3.56 


1 2.05 


46 


0.38 


3.68 


0.90 


1 .453 E 7 


1 8.0 


9.843 


4.807 


0.823 


1 9.53 


1 6.21 


46 


0.52 


3.68 


0.90 


1 .954 E 7 


1 9.6 


9.842 


4.809 


0.823 


28.27 


22.26 
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0.90 
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4.807 
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0.90 


6.538 E 7 


26.5 


9.773 


4.809 


0.819 


122.4 


73.81 
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0.90 


8.743 E 7 


28.7 


9.745 


4.808 


0.817 


176.4 


99.28 
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3.68 


0.89 


1.169 E 8 
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9.700 


4.809 


0.814 


254.5 
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0.89 


1 .547 E 8 
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9.639 


4.808 


0.809 
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0.89 


2.049 E 8 
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4.809 


0.804 


529.0 


240.3 
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0.89 


2.690 E 8 
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9.447 


4.807 


0.796 
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321.5 
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10 


3.69 


0.89 


3.502 E 8 


43.8 



water for given Ra and Pr. More than 330 data points are plotted on figure 10, eacli Nu 
lias been averaged over typically 1 hour and after a relaxation time of a few hours [tne is 
a few tens of seconds) . Obtaining so many data would have required about one year and 
half of continuous operation in water. If the bimodality results from a macroscopic degree 
of freedom, it is likely that the observation of two modes was possible in our experiment 
thanks to the small times scales involved. 
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Let us now present a consequence of this bimodality for room temperature convec- 
tion cells, assuming that bimodality results from the slow-dynamics degree of freedom. 
In heat transfer measurements or visualisation, the possibility to be in a slow transient 
mode-flipping should be considered. The resulting transient-regime uncertainty may be 
a delicate experimental issue for flow characterization studies, and a limiting factor for 
precise measurements, at least without new specific cell design. 



15 



p 


T 


P 


Q 


AT 


ATadiab 


a.AT 


W 


Pr 


Ra 


Nu 


IKg/m-'j 


w 




















51 .29 


5.844 


3.07 


2.01 1 


0.7785 


66 


0.088 


2.80 


2.5 


9.386 E 7 


30.8 


51 .28 


5.845 


3.07 


3.236 


1 .1 08 


66 


0.1 2 


2.80 


2.5 


1 .371 E 8 


34.4 


51 .27 


5.845 


3.07 


5.202 


1 .571 


66 


0.18 


2.80 


2.5 


1 .978 E 8 


38.9 


51 .25 


5.845 


3.07 


8.352 


2.252 


66 


0.25 


2.80 


2.5 


il.OOD E O 


43.6 


51 .22 


5.845 


3.07 


1 3.41 


3.247 


66 


0.37 


2.80 


2.5 


4.1 66 E 8 


48.8 


51.19 


5.847 


3.07 


21.51 


4.687 


66 


0.53 


2.80 


2.5 


f r\'{ o i~ o 

6.01 8 E 8 


54.7 


51.18 


5.847 


3.07 


21.50 


4.705 


66 


0.53 


2.80 


2.5 


6.034 E 8 


54.4 


51.18 


5.846 


3.07 


24.46 


5.1 68 


66 


0.58 


2.80 


2.5 


6.640 E 8 


56.5 


51.15 


5.847 


3.07 


30.80 


6.1 99 


66 


0.69 


2.80 


2.5 


7.956 E 8 


59.6 


51.13 


5.848 


3.07 


34.49 


6.781 


66 


0.76 


2.80 


2.5 


8.68b E 8 


61 .1 


51.12 


5.848 


3.07 


38.77 


7.430 


66 


0.83 


2.80 


2.5 


r" o o r o 

9.522 E 8 


62.9 


51 .08 


5.848 


3.07 


49.05 


8.921 


66 


1 .00 


2.81 


2.5 


1 .140 E 9 


66.6 


51 .05 


5.849 


3.07 


55.30 


9.798 


66 


1 .1 


2.81 


2.5 


T 't AO c r\ 

1 .248 E 9 


68.5 


51 .03 


5.849 


3.07 


61.40 


1 0.61 


66 


1 .2 


2.81 


2.5 


1 .351 E 9 


70.5 


50.98 


5.850 


3.07 


77.09 


1 2.71 


66 


1 .4 


2.81 


2.5 


1 .609 E 9 


74.2 


50.94 


5.851 


3.07 


88.65 


1 4.21 


66 


1.6 


2.81 


2.5 


1 .791 E 9 


76.6 


50.91 


5.851 


3.07 


97.84 


1 5.29 


66 


1 .7 


2.81 


2.5 


1 .924 E 9 


78.8 


50.83 


5.853 


3.07 


1 23.5 


1 8.43 


66 


2.0 


2.81 


2.4 


T TOO C O 

2.298 E 9 


83.0 


50.77 


5.854 


3.07 


1 42.1 


20.58 


66 


2.3 


2.82 


2.4 


2.549 E 9 


85.7 


50.74 


5.855 


3.07 


1 54.9 


22.00 


66 


2.4 


2.82 


2.4 


2.71 7 E 9 


87.7 


50.62 


5.857 


3.07 


1 96.1 


26.58 


66 


2.9 


2.82 


2.4 


3.240 E 9 


92.3 


50.53 


5.859 


3.07 


227.8 


29.98 


66 


3.3 


2.83 


2.4 


3.61 b E 9 


95.4 


50.48 


5.859 


3.07 


247.3 


32.1 


66 


3.5 


2.83 


2.4 


3.859 E 9 


96.9 


50.33 


5.863 


3.07 


310.4 


38.33 


66 


4.1 


2.83 


2.4 


4.519 E 9 


102 


50.18 


5.866 


3.07 


365.0 


43.94 


66 


4.7 


2.84 


2.4 


5.093 E 9 


106 


50.1 3 


5.867 


3.07 


391 .5 


46.30 


66 


4.9 


2.84 


2.4 


c o o o c r» 

5.338 E 9 


1 08 


49.69 


5.875 


3.07 


584.9 


64.1 3 


66 


6.7 


2.85 


2.3 


7.035 E 9 


117 


49.37 


5.883 


3.07 


736.8 


77.60 


66 


8.0 


2.87 


2.3 


8.1 81 E 9 


1 23 


48.98 


5.892 


3.07 


937.2 


94.60 


66 


9.5 


2.88 


2.3 


9.521 E 9 


1 29 


48.52 


5.902 


3.07 


1 1 80 


1 1 5.0 


65 


1 1 


2.90 


2.2 


1 .094 E +1 


1 35 


47.94 


5.91 5 


3.07 


1 502 


141.1 


65 


1 3 


2.92 


2.2 


1 .254 E +1 


1 41 


47.29 


5.931 


3.08 


1 891 


1 72.5 


65 


1 6 


2.95 


2.1 


1 .41 o E +1 


1 46 


51.34 


5.572 


2.72 


1 .239 


0.4834 


70 


0.083 


2.67 


3.4 


1 .1 72 E 8 


33.8 


51 .33 


5.572 


2.72 


2.000 


0.6756 


70 


0.1 2 


2.67 


3.4 


1 .71 5 E 8 


37.7 


51 .32 


5.573 


2.72 


3.21 7 


0.9666 


70 


0.1 7 


2.67 


3.4 


1 coo r" o 

2.533 E 8 


41 .5 


51 .30 


5.573 


2.72 


5.1 69 


1 .351 


70 


0.23 


2.67 


3.4 


3.61 2 E 8 


47.5 


51 .31 


5.573 


2.72 


5.1 71 


1 .359 


70 


0.23 


2.67 


3.4 


o o o r" o 

3.638 E 8 


47.2 


51.31 


5.573 


2.72 


7.900 


1.859 


70 


0.32 


2.67 


3.4 


5.036 E 8 


52.6 


51 .28 


5.573 


2.72 


7.886 


1 .861 


70 


0.32 


2.67 


3.4 


5.037 E 8 


52.4 


51 .29 


5.573 


2.72 


7.854 


1 .858 


70 


0.32 


2.67 


3.4 


5.038 E 8 


52.3 


51 .30 


5.572 


2.72 


9.728 


2.1 74 


70 


0.37 


2.67 


3.4 


5.935 E 8 


55.4 


51 .28 


5.573 


2.72 


9.901 


2.1 99 


70 


0.38 


2.67 


3.4 


5.979 E 8 


55.8 


51 .28 


5.573 


2.72 


1 2.21 


2.599 


70 


0.44 


2.67 


3.4 


7.1 1 8 E 8 


58.2 


51 .26 


5.574 


2.72 


1 2.48 


2.632 


70 


0.45 


2.68 


3.4 


7.1 72 E 8 


58.8 


51 .24 


5.573 


2.72 


1 3.32 


2.781 


70 


0.47 


2.68 


3.4 


7.584 E 8 


59.4 


51 .27 


5.573 


2.72 


1 5.41 


3.079 


70 


0.53 


2.67 


3.4 


8.451 E 8 


62.2 


51 .24 


5.573 


2.72 


1 5.78 


3.1 54 


70 


0.54 


2.68 


3.4 


o.o24 E o 


62.2 


51 .24 


5.573 


2.72 


1 9.46 


3.669 


70 


0.63 


2.68 


3.4 


1 .007 E 9 


66.2 


51 .20 


5.574 


2.72 


21.37 


4.002 


70 


0.68 


2.68 


3.4 


1 .093 E 9 


66.5 


51.19 


5.574 


2.72 


24.52 


4.430 


70 


0.75 


2.68 


3.4 


1 .21 1 E 9 


69.2 


51.1 5 


5.574 


2.72 


30.87 


5.303 


70 


0.90 


2.68 


3.4 


1 .447 E 9 


73.0 


51.16 


5.574 


2.72 


30.87 


5.309 


70 


0.90 


2.68 


3.4 


1 .448 E 9 


72.9 


51.13 


5.575 


2.72 


34.26 


5.779 


70 


0.98 


2.68 


3.4 


1.571 E9 


74.4 


51.1 1 


5.575 


2.72 


39.05 


6.388 


70 


1.1 


2.68 


3.4 


1.737 E 9 


76.9 


51.06 


5.576 


2.72 


49.37 


7.657 


70 


1.3 


2.68 


3.4 


2.070 E 9 


81.5 


51.00 


5.576 


2.72 


62.46 


9.235 


70 


1.6 


2.68 


3.4 


2.479 E 9 


85.9 


50.93 


5.577 


2.72 


79.02 


11.09 


70 


1.9 


2.69 


3.4 


2.951 E 9 


90.9 


50.84 


5.579 


2.72 


100.0 


13.34 


70 


2.2 


2.69 


3.3 


3.507 E 9 


96.2 


50.77 


5.579 


2.72 


123.8 


15.69 


70 


2.6 


2.69 


3.3 


4.092 E 9 


102 


50.75 


5.580 


2.72 


123.7 


15.76 


70 


2.6 


2.69 


3.3 


4.095 E 9 


101 



4.2 Side- Wall conductivity 

A second outcome of this experiment was the finding of a significant and unexpected 
infiuence of the side-wall conductivity on the apparent Nusselt number, as also discussed 
in references [42, 43, 44] : typical side- walls found in the literature are responsible for 
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overestimation of Nu up to more than 20%. We have conducted a specific study, both 
theoretical and experimental for thin wall cells[ll]. 

As discussed in ref. [11], the influence of thin walls can be characterized by a dimen- 
sionless number, called the side-wall number W. This number is the ratio between the 
conductance of the empty cell and that of the fluid at rest. Typical values of W for ref- 
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crcncc experiments range from nearly up to 3.5. Points on figure 11 gather our heat 
transfer measurements (restricted to the lower "modality"). 

The starting point are the same as reference [42] : the wall is in contact with a fluid 
of nearly uniform temperature T through a lateral boundary layer. In contrast with ref- 
erence [42], we do not take the the thickness of the boundary layer as constant (model 2 
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of [[42]]) nor proportional to l/\fRe (model 1 of [[42]]) but proportional to the thermal 
boundary layer of the plates, that is 1/Nu. It results in an exchange length between the 

wall and the bulk proportional to 1/V Nu. Another difference with reference [42] is that 
we consider this exchange length as extending the effective area of the plate, instead of 
considering the calculated flux as a thermal leak. This point of view is clearly confirmed 
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Figure 10: NuRa'^-^^ vs. Ra in the mini-cell. • Pr ^ 0.73 - 0.76 □ Pr ~ 0.82 - 0.84, O 
Pr ~ 0.89-0.90, ■ Pr ~ 0.85-1.1 x Pr ~ 1.6-1.9, a Pr ~ 2.1-2.5, S Pr ~ 2.6-2.7, 
□ Pr ~ 2.9 - 3.4, ♦ Pr ~ 3.6 - 4.3, ■ Pr ~ 5.9 - 6.0, A Pr ~ 11, T Pr ~ 14 - 15, For 
comparison, large cell[8] data are plotted with the symbols : A Pr ~ 0.7, O Pr ~ 1.1. 

by numerical simulation [43]. A Nu correction based on the leak, rather than corrected 
exchange surface, gives a poorer fit on the experimental data (less of 2 decades of Ra are 

well corrected). Our picture result in an asymptotic correction in y/W fully confirmed by 
our compilation 11. 

The closed analytical correction formula we derived has an adjustable parameter ac- 
counting the proportionality between the lateral boundary layer thickness and the plates' 
thermal ones. 

We derived a closed-analytical correction formula for the Nusselt number with one 
adjustable parameter (continuous line of figure 11). Calling Numea and Nucor the measured 
and corrected Nusselt numbers : 

NUmea = Nu,or{l + f{W)) (8) 

with 



and 

^ = 0.8 (10) 

The dependence of this correction with the Nusselt itself has been validated for Rayleigh 
numbers covering 4 decades (and for various values of W) . On the figure, the vertical arrows 
represent the correction magnitude numerically estimated by Verzicco[43]. 
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Figure 11: Nu.Ra~^'^^ vs. the square root of the side-wall number W. The symbols are 
experimental data obtained in the lower mode, for 10® < Ra < 5.10® mh = 2 cm and 20 cm 
cells. The vertical bars represent the magnitude of the side-wall correction estimated from 
numerical simulations by Verzicco[43]. The solid line is the present analytical model with 
one adjustable parameter. 

The present analysis has numerous consequences. First it explains some surprising 
results and discrepancies between published results, as detailed in ref. [11]. Second, the 
correction on Nu being Ra dependent, it changes the Nu{Ra) apparent exponent, at least 
for Ra < 10^°. The examination of several published results indicated that the measured 
exponent -often close to 2/7- are significantly underestimated. Values larger than 0.3 are 
obtained after correction (some experiments claiming a 2/7 exponent are not subject to 
side- wall correction [21]). Thirdly, it should be noted that the side- wall effect mimics a 
Ra dependent exponent, which could be falsely interpreted as an indication of non-power 
law behaviour of convection. Finally, the wall-fluid interaction is likely to introduce a new 
length scale in the convection problem if the wall thickness is non-uniform (flanges, large 
o-ring,...). Such artefact could also cause apparent non-power law behaviour. 

On the other hand, we can consider as very good news that such an important effect 
can be corrected. It could modify the flow itself in such a way that no comparison would 
be possible between the unperturbed and the perturbed case. As shown by Verzicco[43], it 
is indeed the case when the conductance of the wall is too high. Our experimental study, 
presented on fig. 11, shows that the wall number well correlates the various data which 
shows the pertinence of the correction. 

All the data presented in this paper are side-wall corrected. Note that the magnitude 
of the bimodality presented above is not affected by this correction. 



4.3 Prandtl number dependence 

Our experiment has been designed to study the influence of Prandtl number near the 
diffusivity cross-over Pr ~ 1. We find a very small -if any- Prandtl number dependence 
over 1.5 decade[45]. This dependence corresponds to an exponent smaller (in absolute 
value) than 0.03 in a power law picture. This result is compatible with ref. [13] but not 
with ref. [12]. Also, the 2/7 theories[46, 47] predict an exponent —1/7 ~ —0.14 which is 
incompatible with our result. 
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Figure 12: Nu.Ra~^/^ versus Pr with -from bottom to top- Ra = 10^'^^ (light gray), 
Ra = 10^ (dark gray) and Ra = 10^° (black). The value Ra = 10^'^^ corresponds to the 
beginning of the hard turbulence regime for aspect ratio F = 0.5. • present data in the 
lower mode, ■ present data in the upper mode, o © data from ref. [13], F = 0.5 and 1 ( © 
corresponds to Ra = 10^), DBS data from ref. [17] , F = 1 (a : Ra = lO^"^^, B : Ra = 10^). 
The lines are from Grossmann and Lohse[48] model with adjustable parameters tuned by 
these authors on the o © data restricted to aspect ratio F = 1. 

A comparison with the prediction of Grossmann and Lohse[48] (G.L.) can be performed 
using their 5 fitting parameters adjusted to fit the data of ref. [13]. It should be emphasized 
that the Ra and Pr overlap between the data of ref. [13] and ours is limited as shown on 
figure 12. Consequently we are testing a prediction of the G.L. theory : an extrapolation 
on the lower Pr side. We find a good agreement since this prediction falls within the error 
bar of nearly all of data (figure 12). We should mention here another recent test [17] of G.L. 
theory on the higher Pr side. These data are also in good agreement with the prediction 
as shown on figure 12. These three set of data are for two different aspect ratio (F = 0.5 
and 1), but once the wall effect corrected, the influence of the aspect ratio seems weak[42]. 
In particular, all the set of data from the different groups suggest that the transition from 
the low Prandtl regime to the high Prandtl one occurs in the neighbourhood of Pr ~ 1 
and not 0.1 as proposed by Kraichnan[40, 9]. 



5 CONCLUDING REMARKS 

We have measurement Nu{Ra) and Nu{Pr) dependences which are both incompatible 
with the 2/7 and 1/3 theories [40], at least under their present form. The Grossmann and 
Lohse theory [48] can account for our data but the discriminating testing of the 5 fittings 
parameters have to be made on a larger range of Ra and Pr numbers. 

The bimodality effect indicates that the mean flow confinement has a significant in- 
fluence (up to few percents) on the precise Nu{Ra) dependence and this influence should 
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hold for all aspect ratio of order 1. Since confinement effects (multi-modality, poor spa- 
tial homogeneity on boundary layers, ...) are not considered by present theories, their 
predicting power is indeed limited in precision. This darkens the perspective that very 
precise NulRa) measurements can discriminate between competing theories. This conclu- 
sion is reinforced by the boundary conditions influence on the global heat transfer, such 
as side-wall conductivity [11, 42, 43] and hole-burning effects in plates[49, 22, 23]. This 
underlines the importance of alternative approaches to probe the heat transfer mechanism 
and (in)validate theories. It also calls for a new generation of cell design with a specific 
attention dedicated to the mean flow and thermal boundary conditions. For example, the 
influence of the large scale flow on the heat transfer suggests that large aspect ratio cell 
could be required to observe true power law scalings. 
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